Thermodynamics of Lemaitre — Tolman — Bondi Model 



Subenoy Chakraborty*, Nairwita Mazumder', Ritabrata Biswas*. 

1 Department of Mathematics, Jadavpur University, Kolkata-32, India. 
(Dated: June 15, 2010) 

Here we consider our universe as inhomogeneous spherically symmetric Lemaitre ~ Tolman — 
Bondi Model and analyze the thermodynamics of this model of the universe. The trapping horizon 
is calculated and is found to coincide with the apparent horizon. The Einstein field equations are 
shown to be equivalent with the unified first law of thermodynamics. Finally assuming the first law 
of thermodynamics validity of the generalized second law of thermodynamics is examined at the 
apparent horizon for the perfect fluid and at the event horizon for holographic dark energy. 
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INTRODUCTION 



The discovery of Hawking radiation [1] completes the cyclic of identifying black hole (BH) as a thermo- 
dynamical object — the laws of BH physics and thermodynamical laws are equivalent. Since then there 
q 1 is a series of works [2] dealing with thermodynamical studies of the universe as thermodynamical system. 
H ' Considering homogeneous and isotropic FRW model of the universe, most studies deal with validity of the 
generalized second law of thermodynamics (GSLT) starting from the first law when universe is bounded by 
bJ). the apparent horizon [3]. Considering various matter system and different gravity theories, it is generally 
l— "™ '■ found that there is a nice agreement of the thermodynamical laws with apparent horizon as the boundary. 
I \ Also it is found that first law of thermodynamics and (modified) Einstein equations are equivalent at the 
apparent horizon . In contrast, there are few works [4] dealing with thermodynamics of the universe with 
Q\ ' event horizon as the boundary. Due to existence of the event horizon, the matter here is chosen as either 
t-H , quintessence or exotic in nature. Here validity of GSLT put some restrictions either on geometry or on the mat- 
ter itself except when the matter is in the form of holographic dark energy (HDE) [5] , no constraint is necessary . 

(N ; 

\& • In the present work, we consider our universe as in homogeneous Lemaitre — Tolman — Bcmdi[LTB] Model. 
This simple inhomogeneous cosmological model agrees with current supernova and some other data [6]. Also 
very recently Clarkson and Marteens [7] give a justification for inhomogeneous model from the point view of 
perturbation analysis. The apparent horizon and the trapping horizon coincide for the model. We are able 
^ . to show that Einstein field equations and unified first law are equivalent on the apparent horizon. Finally we 
determine the constrains to satisfy the GSLT on the apparent horizon for the perfect fluid and on the event 
• horizon with matter as HDE. 

a- 

_ _ i 

II. BASIC EQUATIONS IN LTB MODEL: 

The metric ansatz for inhomogeneous spherically symmetric LTB space time in a co-moving frame is given 

by 



dS 2 = -di 2 + R ' 2 , dr 2 + R 2 (d8 2 + sin 2 8dcj) 2 ) (1) 
1 + f{r) 

where R = R(r, t) is the (area) radius of the spherical surface and f(r) (> —1) is the curvature scalar (classifies 
the space-time as bounded, marginally bounded and unbounded depending on the range of its values which are 
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respectively /(r) < 0, f(r) = 0, f(r) > 0). Let us suppose that the universe is filled with perfect fluid with 
energy momentum tensor 

Tfiu = {P + p)unU v + V9^v ■ (2) 

where p and p are respectively the matter density and pressure of the fluid and u M is the fluid-four velocity 
of the fluid with normalization u^u v = — 1. By introducing the mass function F(r,t) [8] (related to the mass 
contained within the co-moving radius r)as 

F(r,t) = R(R 2 f(r)) (3) 

the Einstein equations can be written as 

R 2 R' R 2 R V ; 

and the evolution equation for R is 

2RR + R 2 + 8nGpR 2 = f(r) (5) 
The energy momentum conservation relation T£. v = gives 

p + 3H(p + p) = and 

p' = (6) 
where H = | + 2-^ , is the Hubble parameter. 

Further the LTB line element can also be written as 

ds 2 = h ab dx a dx b + R 2 dn 2 (7) 

where 

dtt% = d6 2 + sin 2 8d(j) 2 is the metric on unit two sphere 

and 

h a b = diag I — 1, — — I , (a, b = 0, 1 with x° — t,x x = r) 

V 1 + f(r)J 

is the metric on the 2D hyper surface normal to the 2-sphere. We now introduce two null vectors <9+ and d- 
normal to the 2-sphere (i.e. on the 2D hyper surface) as 



d + = -V2 (ft - Vl+ R ; f{r) d r j and d_ = -y/2 (d t + ^/^ j (8) 

The dynamical apparent horizon which is essentially the marginally trapped surface with vanishing expansion 
is a spherical surface of radius R = Ra satisfying [9] 

h ab = d a Rd b R = 

i.e., 

R A = F{r, t) and R 2 A = 1 + f(r) (9) 
A trapping horizon (Rt) is defined as a hyper surface foliated by marginal spheres and is characterized by [9] 

d+ R T = 



i.e. Rt = y/1 + f(r) (10) 

Thus trapping horizon coincides with the apparent horizon and the result is in agreement with Lemma (I) 
inref.[9]. 



III. THERMODYNAMICS ON THE APPARENT HORIZON: 

We start this section with the following theorems: 

Theorem 1. The unified first law is equivalent to the Einstein field equations at any spherical 
surface. 

Proof:- The unified first law states that [10] 

dE = Aip + WdV (11) 

where, 

R 



is the Misner-sharp mass. 



E = ^ (1 - h ab d a Rd b R) , 
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* = ip a dx a , A = 4irR z , V = -ttR\ areal volume 



i> a = T b a d b R + Wd a R, 
is the energy flux (or momentum density) and 

W = -^Trace{T), 

is the work function or energy density Note that here trace is referred to the two-dimensional space normal to 
the spheres of symmetry. 

For the present LTB model, 

E=^(R 2 -f(r)) (12) 



So 



dE = ^[R{R 2 + 2RR- f{r))dt + {R\R 2 - f(rj) + R(2RR' - f'(r)))dr] (13) 



W = \{p-p) 



V^o = -\{p + P )R , in = \{ P + P)R' (14) 



Therefore, 



Thus 



* = - - (p + p) (Rdt - R'dr) 



AtJj + WdV = AttR 2 pR'dr-pRdt (15) 
Hence equating (13) and (15) according to the unified first law (11) we have 

R 2 + 2RR - f(r) = SirGpR 2 (16) 



4 



and 

(R'(R 2 - f(r)) + R{2RR> - /'(r))) = 8ttGR 2 R' P (17) 

We see that eq. (16) is nothing but the evolution equation (5). Also equation (17) can be written (after some 
simplification) 

£{{#-mR)R} F ^ t) 

8lTGp= WW 

The other equation in cq. (4) for p can be obtained by differentiating equation (3) and using the evolution 
equation (5). Therefore, we write 

Unified First Law of Thermodynamics Einstein Equations. 

One may note that in rcf [10], Cai et al stated that unified first law is an identity concerning the (0, 0) 
component of the Einstein equation. But here we have shown that unified first law is more general — it is 
equivalent to the Einstein field equations at any spherical surface of symmetry. 

Theorem 77. The validity of Clausius relation at any spherical surface of symmetry depends 
on the choice of the tangent vector. 

Proof : The Clausius relation states that 

z)= £g {dA z) (18) 

where k is the surface gravity and z is any vector tangential to the spherical surface of symmetry Let us choose 

z = z + d + + z~d- 



where z\, Z2 (i.e., z + , z~) are constant parameters. By definition, the surface gravity is given by 

k = —^=d a (V^h h ab d b R) , 

2y/—h V / 

which for the present model has the expression 



(19) 



VTT7 

2R' 



R" 



R' R + RR' 



v/T+7 



(20) 



Now, {Aip, z) = {dA, z) gives on simplification, 

zi = yr+7 

z 2 R 



2ttR jp + p) - | 
2irR (p + p) + g 



(21) 



This clearly shows that the ratio ^ will have different values at different spherical surface. Therefore, validity 
of the Clausius relation depends on the choice of the approximate tangent vector and is in agreement with 
lemma III of ref [9] . 

Now assuming the Clausius relation at the (apparent or event) horizon we examine the validity of the GSLT. 
Let Sh and Si be the entropy of the horizon and the matter bounded by the horizon. Then GSLT states that 



^ {S H + St) > 



(22) 



Since we are considering equilibrium thermodynamics so temperature of the inside matter distribution can be 
considered as that of the horizon {T H ). 
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Now to find the variation of entropy at the horizon we start with unified first law 

dE = + WdV 

using equation (14) and (15) we write 

dE = -Anr 2 (p + p) Rdt + 4ttR 2 pR'dR (23) 
Hence from the Clausius relation at a horizon 

T H dS H =SQ = -dE = 4ttR 2 h (p + p) R H dt 

i.e., 

d_S» = 4,R 2 H R H (p + p) dt 
dt Tfj 

where Sh and Tu are the entropy and temperature at the horizon. 

In order to determine the time variation of the matter entropy. We start with the Gibb's equation [11] 

T H dSj = dE t + pdV (25) 

where V — ^itR^, Ej = pV and (p, p) are the matter density and pressure of the fluid bounded by the horizon. 
Now using the energy conservation equation (6) we have from Gibb's equation (after some simplification) 

T H ^ = 4nR% (P + P)^- ^R'hH (p + p) (26) 

Hence combining (24) and (26) we have, 

m d ,„ \ 4ttR 3 h , . i 4R H R' H ) . 

T H -(S B + Sl ) = ^.{p + P)[-^-^ (27) 

Now we shall examine the validity of GSLT both at the apparent horizon and at the event horizon. 
Case-/ : Universe filled with perfect fluid and bounded by the apparent horizon. 

The apparent horizon for the LTB model is characterized by 

R 2 = 1 + f(r) 

or equivalently, 

R = F. 

So eq. (27) now becomes 



TA dt {SA + Sl) - 3F'(r, t) 



F' f(r) 



4— 



F 2{l + /(r)} 



(28) 



IV. CASE-// : UNIVERSE FILLED WITH HDE AND BOUNDED BY THE EVENT HORIZON 

As geometrically event horizon (Re) can not be evaluated for LTB model so we try to evaluate Re from 
physical consideration. It is HDE in which energy density (pd) can be written as [12] 



E 



G 



where, c is any free dimensionless parameter estimated by observational data [13]. However , in the present 
work we have taken c to be arbitrary. Now from the energy conservation relation (6) we obtain 



Re = -HR E (l + u D ) 



where ujd 



^ is the equation of state parameter for the HDE. Then from cq (27) 



~dt {Se + Sl) ~ W E 



3 R 'e 
Re 



1 + 0J D 



(30) 



(31) 



Now we shall analyze the above results for the validity of GSLT (i.e., inequality (22)). 
Apparent Horizon 
(a) Quintessence Era (p + p > 0) 

In an expanding universe both R and R' are positive. So from the expression of p (in eq (4)) we see that 
F 1 must be positive. But F is negative or positive depending on whether p is positive or negative. Also from 
equation (3) F should be positive. Hence from equation (27) GSLT will be valid if 

«i» < o 



or 



m>0 and ^(m^ 



> 



(b) Phantom Era (p + p < 0) 

In this case GLST will be valid if 



f(r)>0 and |-(ln-^)<0 



Therefore, validity of GSLT at the apparent horizon depends on the arbitrary integration functions appear 
in the Einstein field equations. Finally, one may note that here we have not used any explicit expression for 
entropy and temperature at the horizon. 

Event Horizon 

From eq.(31) we see that validity of GSLT does not depend explicitly on whether the HDE satisfies the weak 
energy condition or not. Essentially GSLT will be satisfied at the event horizon if both Re and jjjj^j^j are 
simultaneously increasing (or decreasing) function of r. 

For future work, it will be interesting to examine whether Bekenstein entropy and Hawking temperature 
formulae hold at the apparent horizon for the present inhomogeneous LTB model. Also it will be nice to make 
an attempt for determining event horizon in this model. 
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